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"^An  analysts  is  presented  of  the  combustion  of  graphite  in  a 

high  speed  stream  of  dissociated  air,.  Many  features  of  the  analysis 

/ 

are  quite  gene-al  and  may  be  applied  to  the  oxidation  ot  different 
materials  in  arbitrary  chemically  reactive  environments.  However, 
because  of  the  current  interest  in  the  hyp  rsonic  leading  edge 
problem,  numerical  results  are  presented  here  which  are  directly 
applicable  to  surface  oxidation  at  the  leading  edge  region  of  fins 
and  wings,  and  the  nose  cap  of  axially-symmetric  hypersonic 
vehicles. 

-^The  reaction  rate  controlled- regime  and  transition-regime 
are  ftret  considered  at  length,  and  it  is  shown  how  the  grade  of 
graphite  and  its  specific  chemical  properties  influence  the  over-all 
uxioaiion  rate. 

It  is  then  shown  how  the  governing  equations  of  change  may 

be  reduced  to  a  coupled  set  of  non-linear  differential  equations  of 

the  fifteenth- order  with  variable  coefficients  and  split  boundary 

conditions.  These  differential  equations  are  then  utilized  in  treating 

the  laminar,  compressible,  multicomponent,  chemically  reacting 

boundary  layer  in  the  diffusion  controlled  regime,  and  solutions 

are  obtained  for  both  the  equilibrium  and  frozen  gas  flow  chemical 

constraints.  Wumeeienl  rasulta  are  obtained  on  an  IBM  704  Ut^'^ral 

computer,  amt  correlated  results  are  obtained  for  the  heat  transfer- 
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the-masa  transfer  rate  and  the  skirt  friction  coefficient  for  the 
complete  range  of  suborbital  hypersonic  flight  conditions  of  interest 
In  addttton;  ~irt  older  to  establish  a  better  understanding  ef 
-the-complex  uhysit ‘^.hemical  processes  which  orriiPt  many  details 
of  the  sti'uciiire  of  the  boundary  layer,  including  the  variation  of 
macroscopic  gas  velocity,  gas  temperature,  chemical  composition 


and  chemical  source  terms  through  the  boundary  layer,  are 


presented 


I.  INTRODUCTION 


In  considering  the  design  of  hypersonic  lifting  vehicles,  special 
attention  must  be  given  to  the  leading  edge  surfaces  which  are  exposed 
tij  sustainefi  aerodynamiv.  heating,  and  hence  must  function  for  long  time 
periods  at  leading  edge  temperatures  in  the  vicinity  of 
.  (Refs.  !-6)'. 

A  class  of  superior  carbonaceous  materials  known  as  graphite 
immediately  suggests  itself  becauee  this  form  of  carbon  is  a  refractory 
material  having  high  thermal  shock  resistance,  good  high  temperature 
strength,  excellent  machinabillty,  high  thermal  conductivity,  a  high  subli¬ 
mation  temperature  and  a  relatively  low  oxidation  rate  (see  Table  1). 

The  type  of  graphite  which  is  in  current  use  in  industry  is  usually 
manufactured  from  carbon  base  materials,  rather  then  mined  as  the  natural 
substance,  and  hence  is  commonly  called  "artificial  graphite".  The  latter 
is  superior  to  either  natural  graphite  or  carbon,  both  of  which  have  relatively 
low  mechanical  strength.  It  is  noted  that  manufactured  graphite  is  not  one 
specific  material,  but  a  family  of  materials  which  are  all  essentially  pure 
carbon.  They  differ  from  each  in  other  in  regard  to  the  orientation  of  the 
crystallites,  grain  s'zip,  the  size  and  number  of  pore  spaces,  the  degree 
of  graphitization,  and  the  level  cf  impurities,  which  therefore  lead  to  certain 
differences  in  the  physical  and  chemical  properties.  Thus,  the  wide  variation 
found  in  the  properties  of  graphite  can  be  attributed  to  the  choice  of  starting 
materials,  and  to  the  degree  of  control  during  the  manufacturing  process.  In 
reference  7  the  reader  will  find  a  concise  review  of  the  properties  and 
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applicatioaa  of  different  grades  of  manufactured  graphite. 

Chemical  reactions  between  carbon,  coal,  graphite  and  oxygen 
have  been  studied  intensively  for  over  one  hundred  years,  and  attention 
has  been  given  to  the  reaction  rate  controlled,  transition,  and  diffusion 
controlled  oxidation  regimes.  Consequently,  a  voluminous  literature  exists, 
and  excellent  reviews  on  the  mechanism  of  carbon  oxidation  have  been 
written  byGolcvina  (Ref.  «),  Frank  -  Kamenetskii  (Ref.  9) ,  Arthur  (Ref.  10), 
Townsend  (Ref.  U).  Strickland  -  Constable  (Ref.  1^) ,  von  Loon  and  Smeets 
(Ref.  1  i) .  Gerstein  and  Coffin  (Ref.  1-1).  Khitrin  (Ref.  IS),  and  Blakeley 
(Ref.  lo).  However,  very  little  of  this  previous  work  appli-ts  to  the  en¬ 
vironmental  conditions  encountered  during  hypersonic  flig.-t,  .Specifif  , 
for  the  regime  of  greatest  interest  to  the  glide  vehicle  designer,  little 
information  is  available  otlier  than  the  theoretical  work  of  Scala  (Refs.  17, 
la),  Lees  (Ref.  19),  Oennison  and  Dooley  (Kef.  ZO),  and  Moore  and 
7.Iotnick  (Ref.  dl). 

None  of  the  previous  work,  either  theoretical  or  experimental,  considers 
the  problem  of  determining  systematically  the  relationship  between  mass  loss, 
heat  transfer,  and  viscous  skin  friction,  as  a  function  of  the  significant 
environmental  parameters,  such  as  the  flight  speed,  the  ambieu.  pressure, 
the  surface  temperature  and  the  .ncdel  geometry. 

In  the  study  presented  here,  the  hypersonic  ablation  of  graphite  is 
considered,  and  the  heat  transfer  and  mass  transfer  processes,  and  viscous 
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drag  effect*;:  which  are  experienced  by  hypersonic  vehicles  flying  in  the 
"earth's  atmosphere,  are  analysed  in  detail.  Although  the  analysis  <ieveloped 
here  is  quite  general,  because  of  great  current  interest,  numerical  results 
have  been  obtained  which  are  applicable  to  the  (leading  edges  of  fins  and 
wings,  and  at  the  forward  stagnation  point  of  axially-symmetric  vehicles. 

Upon  inr-uducing  available  experimental  data  on  heterogeneous 
reaction  kinetics  (Refa'f,  ^2-26)  ,the  mass  transfer  and  heat  transfer  rates 
will  be  determined  at  low  surface  temperatures,  fnr  the  reaction  rate  con¬ 
trolled  regime.  At  higher  surface  temperatures,  the  transition  regime 
behavior  will  be  determined  utilizing  the  results  obtained  in  the  rate  con¬ 
trolled  and  ditfusion  controlled  regimes,  by  applying  the  concept  of  resistances 
to  mass  transfer  in  series.  Since  the  heat  transfer  to  the  surface  depends 
on  the  ratio  of  carbon  monoxide  to  carbon  dioxide  at  the  surface,  the  recent 
data  of  Arthur  (Ref.  27)  and  Bonnetain  (Ref.  28)  are  also  introduced. 

At  still  higher  surface  temperatures,  in  the  diffusion  controlled 
regime,  exact  solutions  will  be  obtained  for  the  laminar  flow  of  a  compressible 
multicomponent  chemically  reacting  gas  over  a  reacting  solid.  It  will  be 
assumed  that  dissociated  air  produced  by  the  upstream  shock  waves  can  be 
treated  as  a  four  component  gas,  consisting  of  oxygen  and  nitrogen  atoms, 
and  oxygen  and  nitrogen  molecules.  Since  the  primary  combustion  products 
include  carbon  monoxide  and  carbon  dioxide,  the  total  number  of  gaseous 
species  considered  is  six.  Therefore,  the  analysis  requires  the  solution  of 
a  coupled  set  of  non-linear  partial  differential  equations,  (including  the 
conservation  of  mass,  momentum,  energy  and  chemical  species)  v.  ^ch  is 
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of  the  .-•I  'l'  r,  ri.ivins;  liouiiUAry  coniUtioiisi  anti  variable 

transport  and  thermodynamic  coefficient*. 

In  the  work  presented  here,  as  in  earlier  atudie*  (Refs.  29.30), 
the  transport  properties  of  the  individual  atomic  and  molecular  species 
Will  be  calculated  utilizing  the  rigid  sphere  and  Leonard- lone ■<  moiel, 
respectively.  The  thermodyiiamic  properties  of  the  pure  species  will  be 
determined  using  the  formulae  of  statistical  mechanics.  The  transport 
and  thermodynamic  properties  of  the  gaseous  mixture  will  be  evaluated 
during  the  solution  of  the  problem,  in  terms  of  the  gas  composition,  the 
pressure  and  the  temperature. 

Since  the  homogeneous  rates  of  reaction  of  the  various  species 
present  fn  the  high  temperature  gas  stream  are  not  yet  known  precisely, 
calculations  will  be  performed  for  the  two  limiting  cases  of  "froBen"flow 
(infinitesimally  slow  cas  phase  reaction  rates),  and  "local  equilibrium." 
flow  (infinitely  fast  gas  phase  reaction  rates),  which  bracket  the  actual 
situation.  It  will  be  shown  that,  as  in  the  case  of  hypersonic  stagnation 
point  heat  transfer  (Refs.  30  to  33i .  when  the  gat  is  in  local  equilibrium 
at  the  surface,  for  arbitrary  hypersonic  free  stream  conditiona,  both  the 
overall  rate  of  rras*  tranafer  and  the  net  heat  transfer  rate  are  virtually 
independent  of  the  rates  of  gas  phase  reaction  in  tha  diffusion  controlled 
regime.  This  precludes  the  necessity  of  having  an  exact  knowledge  of 
gas  phase  kinetics. 

In  order  to  esUblish  a  better  understanding  of  the  physicochemical 
processes,  many  details  of  the  structure  of  the  multicomponent  boundary 
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layer  will  be' presented,  incltidinf;  the  variation  of  velocity,  temperature, 
and  gas  composition  as  a  function  of  distance  from  the  surface.  The  zones 
in  wnioh  chemical  reactions  occur  in  the  gas  phase  and  the  magnitude  of 
tile  various  chemical  source  distributions  w.  will  be  shown  in  detail. 
Calculations  will  also  be  made  to  determine  the  fraction  of  the  heat  trans¬ 
ported  to  the  reacting  surface  by  the  various  fluid  dynamic  and  molecular 
interaction  processes.  Finally,  the  magnitude  of  the  viscous  shear  stress 
will  be  evaluated  and  correlated  in  the  form  of  a  skin  friction  coefficient. 

Utilizing  the  graphical  results,  and  the  correlation  formulas  which 
will  be  presented  here,  one  mav  predict  the  heat  conducted  into  the  solid, 
the  mass  loss  from  tne  leading  edge  and  the  skin  friction  coefficient,  for 
a  wide  range  of  hypersonic  flight  conditions,  i.  e.  Mach  numbers  in  a 
range  from  10  to  24,  an  altitude  range  from  10,000  ft,  to  240,000  ft.  ,  and  . 
surface  temperatures  from  the  threshold  range  through  6000® R,  for  vehicles 
of  arbitrary  nose  radius  or  wing  leading  edge  radius^  angle  of  yaw. 
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II.  .SYMBOLS 


a,b 

Alt. 


C. 

1 


JJ.. 

ij 


coefficients  in  e<]iii1ibrium  constant; 
altitude 


skin  friction  coefticient 


1/2 


/?  U 

j  w  ~e 


mass  fraction  of  species  i 


specific  heat  at  constant  pressure  of  species  i 


7*  CjCp.  ,  frozen  specific  heat  of  the  mixture; 


specific  heat  of  the  solid 


binary  diffusion  coefficient 


self  diffusion  cc'efficient 


multicomponent  diffusion  coefficient 


activation  energy 

similarity  stream  function 

u/u^,  dimensionless  velocity 

static  enthalpy  of  species  i,  including  chemical 
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A  standard  heat  of  formation  of  specie*  i  evaluated  at 


s 


^  C^h^,  static  enthalpy  of  mixture 


A  heat  of  vaporization 


Ji 


diffusion  flux  of  species  i 


mechanical  equivalent  of  heat 


specific  reaction  rate 


effective  collision  frequency 


K 


frozen  thermal  conductivity  of  the  mixture 


Si 


Li 


equilibrium  constant  of  specie*  i 


K 


'•  frozen  Lewis  number 


,  frozen  thermal  Lewis  number 


/’w^'w  *  inte>T>ha*e  mas*,  transfer 


Mj 


molecular  weight  of  species  i 
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mean  molecular  weight  of  the  mixture 


number  of  chemical  species 
oMer  of  the  reaction 

number  of  moles  of  species  i  per  unit  volume 


heat  of  reaction 

energy  transfer  function  =  KVT  Vjh  j 


heat  transfer  rate 

effective  heat  ot  ablation 


universal  gas  constant 
nose  radius  of  body 


Reynolds  number  based  on  x 


solid  state,  condensed  phase 


T 


temperature 


u 


V 


V. 

1 


w. 

1 


y.  r 

y 

A 

€ 

Ti 
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X  component  of  velocity 
y  component  of  velocity 
macroscopic  stream  velocity 
absolute  velocity  of  species  i 
diffusion  velocity  of  species  i 

flight  speed 

chemical  source  term,  mass  rate  of  production  of 

species  i  by  chemical  reaction  per  unit  volume  per 
unit  time 

mole  fraction  of  species  i 
body  oriented  coordinate  system 
linear  rate  of  surface  recession 
yaw  angle 

emissivity,  depth  of  potential  well 
similarity  variables 

dimensionless  temperature 
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characteriatic  vibrational  temperature 
viscosity  coefficient  of  species  i 
viscosity  of  mixture 
density 

Stefan- Bolt^cmann  constant 

viscous  shear  stress 

collision  diameter 

stream  function 

treated  as  if  the  gas  is  dissociated  air 

condensed  phase 

calorimeter 

outer  edge  of  boundary  layer 
equilibrium 

gaseous  species 
i  th  species 
radiation 
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stagnation  point 


vaporization 


wall,  interface 


upstream  of  shock,  edge  of  boundary  layer 


denotes  differentiation  with  respect  to 


m.  THE  REACTION  RATE  CONTROLLED  REGIME 
4)  Mechanism  of  Surface  Degradation 

In  an  investigation  of  the  behavior  of  graphite  in  dissociated  air,  one 
requires  data  on  the  nature  and  extent  of  the  chemical  reactions  between,  carbon 
and  the  primary  products  of  dissociated  air,  including  atomic  and  molecular 
oxygen  and  nitrogen. 

The  recent  work  of  Stieber  (Ref.  i4)  indicates  that  nitrogen  molecules 
can  be  considered  chemically  inert  on  carbon  surfaces  at  temperatures  as  high 
as  5400®  R.  An  early  study  by  Strutt  (Ref.  35)  indicates  that  active  nitrogen 
(primarily  atomic  nitrogen  in  the  ground  state)  does  not  react  with  carbon  at 
room  temperature.  In  addition,  Zinman  (Ref.  36|  who  studied  the  interaction 
between  atomic  nitrogen  a.id  carbon  at  800®  C,  did  not  detect  measurable 
amounts  of  either  cyanogen  or  paracyanogen.  Consequently,  one  may  conclude 
that  at  surface  temperatures  up  to  approximately  2000®  R,  no  permanent 
carbon-nitrogen  compounds  are  formed  at  a  carbon  surface.  While  experimental 
data  is  lacking  at  higher  surface  temperatures,  it  will  be  assumed  that  molecular 
nitrogen  is  chemically  inert,  and  that  atomic  nitrogen  undergoes  heterogeneous 
recombination  at  a  graphitic  surface.  Thus,  if  atomic  nitrogen  diffuses  to  the 
surface  without  undergoing  gas  phase  recombination,  then.the  graphite  acts  to 
catalyse  the  recombination  of  the  atoms  at  the  surface.  In  this  study,  therefore, 
surface  degradation  will  be  assumed  to  be  primarily  a  consequence  of  a  surface 
oxidation  process.  That  is,  chemical  reactions  between  nitrogen  and  carbon,  and 
mechanical  effects  such  as  spalliug,wtIX  not  be  included  in  the  theoretical  calculation 
of  the  ablation  rate,  during  hypersonic  flight. 
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Th*  carbon-oxygM  reaction  hae  been  studiet*  extenaively  for  over 
one  hundred  year*  and  hence,  fortunatelyt  although  the  mechaniam  ia  atill 
not  completely  underatood,  aufiicient  experimental  data  exiata  upon  which 
rcaaonable  theoretical  calculation*  may  be  baaed.  The  fact  that  graphite  burna 
to  form  a  mixture  of  carbon  monoxide  and  carbon  dioxide  hae  h*en  diacuaaed 
in  the  literature,  although  there  ia  aome  diaagreement  aa  to  the  aequence  of 
the  atepa  in  the  chemical  reaction*.  At  thia  time,  aa  wi7'  Sr  diacuaaed,  there 
ia  alao  aome  diaagreement  aa  to  whether  the  overall  reaction  ia  firat  order 
with  reapect  to  the  concentration  of  oxygen  or  of  fractional  order.  It  ia 
intereating  to  note,  however,  that  the  wideat  difference*  in  the  oxidation 
behavior  of  the  varioua  gradea  of  graphite  are  found  at  the  lowest  anrface 
temperature*  where  the  proceaa  ia  rate  controlled,  and  theae  individual 
difference*  tend  to  diaappear  aa  the  aurface  temperature  rise*. 

Since  the  reaction  between  carbon  and  oxygen  produce*  both  CO  and 
C02>  theae  product*  can  be  the  reault  of  either  parallel  or  conaecutive  reactions. 
A  number  of  different  mechanisma  are  possible  and  theae  include  the  following: 

1.  The  formation  of  both  CO  and  CO2  in  a  surface  reaction  between 
C(a)  and  O2  or  O. 

2.  The  formation  of  CO2  from  C(s)  and  O2  or  O  at  the  surface,  followed 
by  the  dissociation  of  CO2  to  CO,  and  O  in  the  gas  phase,  or  reduction  of 
CO2  to  CO  at  the  surface. 

3.  The  formation  of  CO  at  the  surface  from  C(*)  and  O2  or  O,  or 
C(*)  and  CO2,  the  CO  being  oxidised  to  CO2  in  the  ga*  phase. 
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Measurements  c»€  the  gas  composition  in  the  vicinity  of  an  oxidizing 
carbon  surface  have  been  made  by  a  number  of  different  investigators,  including 
Arthur  (Ref.  27),  Bonnetain  (Ref.  28)  and  Snow  et  al  (Ref.  j7).  They  have 
verified  that  both  species  are  detected  adjacent  to  the  surface,  even  nt  low 
surface  temperatures,  and  that  the  ratio  of  the  mass  fraction  of  CO  to  CO2 
at  the  surface  r’^e-;  rapidly  with  increases  in  surface  temperature.  Their 
data  can  be  represented  by  an  Arrhenius  equation, 

(CcoAcOj)^  ■  !»-='«’■.  (1) 

and  are  shown  in  Fig.  '3. 

It  is  noted  that  if  the  gas  at  the  reacting  surface  had  sufficient 
time  to  achieve  thermochemical  equilibrium  during  the  low  temperature 
oxidation  process,  then  the  composition  of  the  gas  could  be  determined 
from  the  equilibrium  constant  for  the  reaction. 


C(8)  +  COp  2CO 


(2) 


Since  the  equilibrium  constant  for  this  reaction  may  be  written  in  the 
form  (see  Table  II), 


(3) 


it  is  not  surprising  tiiat  alien  the  boundary  la, or  aolutions,  which  are 
based  on  the  assumption  of  local  thermochemical  equilibrium  at  the  surface, 
are  compared  with  the  experimental  data  of  .Arthur,  Bonnetain,  and  Snow 
et  al,  see  Fig.  15.  the  same  general  trend  -s  exhibited.  One  may  conclude 
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that  either  the  gaseous  CO-  CO^  system  is  actually  never  removed  very  far 
from  an  equilibrium  state  during  oxidation,  or  that  the  experimental  technique 
utilized  by  the  investigators  produces  a  shift  toward  the  equilibrium  composition. 

Although  these  data  cannot  be  utilized  to  rationalize  the  presence  of 
ary  of  the  three  suggested  mechanisms,  it  is  fortunate  that  the  mass  transfer 
and  heat  transfer  e  surface  can  be  predicted  reasonably  well  without 
specifyinp,  tne  specific  oxidation  mechanism, 
b)  Reaction  Rate  Data 

It  is  commonly  accepted  that  the  manner  in  which  the  oxidation  of 
graphite  proceeds  depends  on  the  type  of  graphite,  the  environmental  conditions 
(e.  g.  the  pressure,  temperature,  velocity  and  composition  of  the  stream),  the 
surface  temperature,  and  at  high  surface  temperatures,  on  the  geometry  of 
the  model.  The  temperature  at  which  a  measuralie  mass  loss  first  occurs, 
is  called  the  threshold  temperature  (1000-1800*^ R)  and  is  not  a  constant  but 
depends  or.  the  partial  pressure  uf  the  reacting  gas  at  the  surface. 

At  luw  surface  temperatures,  the  nttass  loss  increases  rapidly  with 
surface  temperature,  and  the  ablation  rate  is  limited  by  the  speed  of  the 
chemical  processes,  including  adsorption,  reaction  and  desorption. 

At  somewnat  higher  temperatures  (1400-3200°  R),  the  speed  of  the 
chemical  processes  is  comparable  to  the  rate  at  which  fresh  reactant  is  brought 
to  the  surface  and  the  products  of  reaction  are  removed  by  convection  and 
diffusion.  Therefore,  the  overall  process  is  in  a  transition  regime,  where 
the  speed  of  the  overall  oxidation  process  is  limited  by  the  presence  of  two 
resisiaiiucs  lu  series,  one  chemical,  and  the  second  gas  dynamic. 


< 
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At  temperatures  above  3200°  R,  the  chemical  oxidatioa  processes 
are  overshadov'ed  by  the  dynamic  processes.  la  this  diffusioa  controlled 
regime,  the  mass  loss  is  relatively  insensitive  to  the  surface  temperature. 
This  result  has  been  found  experimentally  for  subsonic  flour  (Ref.  37)  and 
will  be  demonstrated  here  for  hypersonic  flow. 

Eventually,  when  the  surface  temperature  is  sufficiently  high  (5500- 
8000°  R),  the  sublimation  rate  of  carbon  atoms  and  molecules  can  exceed  the 
surface  oxidation  rate,  and  these  species  are  then  present  in  the  gas  phase. 

Different  investigators  have  studied  one  or  more  of  these  oxidation 
regimes.  For  example,  the  influence  of  environment  upon  the  combintion 
rate  of  carbon  has  been  studied  by  Hottel  et  al  (Refs.  22,  39,  40),  Chukhanov 
and  Grozdovskii  (Ref.  38),  Gulbransen  (Ref.  42),  and  more 

recently  by  Kuchta,  Kant  and  Darrion  (Ref.  43). 

The  effect  of  the  nature  of  the  carbon  on  the  oxidatioa  rate  has  also 
been  subject  to  investigation.  For  example,  Riley  (Ref.  44),  and  Smith  and 
Polley  (Ref.  45)  have  studied  the  effect  of  varying  degrees  sf  crystallinity  of 
the  carbon.  Winslow  et  al  (  Refs.  46,  47)  and  Akamatsu  et  al  (Ref.  4S), 
investigated  the  relationship  between  the  starting  materials  and  the  degree 
of  graphitization.  Wicke  and  Hedden  (Refs.  49,  S<t>  have  postuI..ted  that  for 
porous  types  of  carbon,  the  diffusion  of  oxygen  into  the  pores  can  be  the  rate 
controlling  step  in  a  transition  regime  between  the  rate  controlled  and  diffusion 
controlled  oxidation  regimes;  hence,  the  apparent  activation  energy  is  found 
to  be  half  the  true  activation  energy.  However,  Blakeley  (Ref.  16),  who 
investigated  natural  and  artificial  graphite  under  a  variety  of  conditions,  feels 
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that  .1  pore  Uifiusion  mechaiiism  is  not  required  to  explain  his  experimental 
cl  it  t 

It  iias  also  beer,  shown,  e.  Arthur  (Refs.  51,  5ii)  ,  that  impurities  . 

;n  ttie  solid  phase,  such  as  sodium  c.xrbonate  and  zinc  chloride,  will  augment  : 
the  rate  of  oxidation  in  the  reaction  rate  controlled  regime.  These  experiments  . 

'.ndicate  that  at  relatively  low  temperatures,  impurities  tend  to  t’.'it  as  favorable 
sites,  ur  catalytic  agents,  which  promote  the  rate  of  reaction.  Although  little  . 
positive  experimental  data  exists,  one  may  anticipate  that  trace  amounts  of 
impurities  in  the  gas  pliase  will  also  influence  the  oxidation  rate  in  the  rate 

•  fiV  * 

controlled  regime. 

With  regard  to  the  dependence  of  the  rate  of  oxidation  upon  the  pressure  . 
of  oxygen  in  the  stream,  there  is  some  experimental  evidence  that  the  reaction 
rate  is  first  order  with  respect  to  oxygen  pressure  (Refs.  22,  49,  53),  However, 
Frar.k-Kamenetskii  (Ref.  9)  has  re-examined  the  experimental  data  of  Parker 
and  Hottel  (Ref.  2.1)  and  has  shown  that  the  reaction  rate  can  be  interpreted  to 
be  of  fractional  order.  Further,  Semechkova  and  Frank-Kamenetskii  (Ref.  54) 
have  shown  that  the  rate  of  reaction  between  carbon  and  carbon  dioxide  in  the 
purely  kinetic  regime  is  lower  than  first  order,  while  Klibanova  and  Frank-Kamenetskii 
(Ref.  55)  have  established  that  the  reaction  between  carbon  and  oxygen  is  not 
first  order,  but  fractional,  lying  between  1/3  and  1/2. 

Vulis  (Ref.  26}  tabulated  a  large  amount  of  data  on  the  kinetics  of  the 
reaction  between  carbon  and  oxygen  or  carbon  dioxide.  Upon  applying  an 
Arrhenius  formula  to  the  data,  he  found  that  the  activation  energy  £  varied 
between  limits  of  8  and  37  Kcal/mble  for  the  carbon-oxygen  reaction,  and  was 
approximately  2.  2  times  larger  for  the  carbon-dioxide  reaction.  Vulis  also 
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found  that  the  logarithm  of  the  specific  reaction  rate  was  a  linear  function 
of  t  ie  activation  enerf;y,and  hence, waa  ledtothe  conclusion  that  the  only 
experimentally  determinable  characteristic  of  a  given  variety  of  carbon  is 
the  activation  energy  of  either  o,  these  reactions.  However,  Vulis'  treatment 
of  the  data  rests  on  the  assumption  that  the  true  chemical  kinetics  at  the 
surface  follow  a  f  .st  order  reaction,  and  it  has  been  pointed  oit  tiy  Vrank- 
Kamenetakii  (Ref.  'll  that  this  assumption  lacks  theoretical  or  experimental 
substantiation. 

Examination  of  a  large  mass  d  experimental  data  (e.  g.  ,  Refs. 

9.  10,  11.  U.  Ui  H.  15.  16,  ^l,  Zi,  Z5,  Z(>,  17,  is,  i9,  40,  41,  Al,  etc., 
indicates  tnat  in  the  reaction  rate  controlled  regime,  the  oxidation  process 
follows  a  rate  law  which  may  be  written  in  the  form: 


■n  =  K(P^  ) 

.  O,  w 

React.  2 


(4) 


where  F*q  is  the  partial  pressure  of  the  element  oxygen  near  the  surface. 


n  is  the  order  of  the  reaction,  and  k  is  the  specific  reaction  rate. 

As  noted,  there  is  some  uncertainty  about  the  precise  value  of  n,  as 
values  in  a  range 


0  n  -S'  l.O  (5) 

have  been  reported  in  the  literature  for  different  oxidation  regimes.  Also, 
as  discussed,  the  specific  reaction  rate  k  is  an  exponentially  increasing 
function  of  temperature  whose  precise  magnitude  is  directly  related  to  the 
type  of  graphite  and  its  treatment  during  manufacture.  Customarily,  *he 

\ 
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reaction  rate  constant  is  written  in  the  standard  Arrhenius  form. 


k  =  k«-E//tT 


(6) 


where  the  p re- exponential  factor  can  vary  over  several  orders  of  maRnitude, 
and  the  activation  energy  has  been  reported  to  fall  within  the  limits 

P  5  E  60  K  cal. /mole  (7) 

For  example.  Gulbransen  and  Andrew  (Ref.  42).  and  Blyholder  and  Eyring 
(Ref.  24)  report  their  data  in  the  form  of  Eq.  (6).  Parker  and  Hottel  (Ref.  22) 
utilized  the  form; 


k  =  k  -E/<«T 

o 

Vulis,  (Ref.  26),  obtained  the  empirical  formula 

k  =  k  x-1^  E(T-T*) /^TT* 
o 

and  Frank-Kamenetskii  (Ref.  9)  has  suggested, 


o 


(8) 


(9) 


(10) 


where  the  associated  rate  data  for  eqs.  (6),  (8),  (9).  and  (10)  appear  in 
Table  III.  These  data  also  appear  in  Fig.  16,  where  it  is  seen  that  the  data 
Parker  and  Hottel,  and  Vulis  are  based  on  a  first  order  reaction,  but  the 
data  of  Gulbransen  and  Andrew,  Frank-Kamenetskii,  and  Blyholder  and 
Eyring,  .have  been  tpVeo  as  following  a  on»-half  order  reaction.  This  figure 
shows  clearly  that  one  can  expect  different  grades  of  graphite  to  exhibit  large 
differences  in  oxidation  behavior  in  the  rate  controlled  regime,  i.  e.  , 
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1400°R  S  ■S  3200“  R.  However,  fortunately  at  higher  surface 

temperatures,  the  process  becomes  diffusion  controlled  and  these  large 
differences  will  no  longer  appear. 

Although  our  subsequent  results  will  usually  be  presented  in  a 
general  form,  wherever  representative  calculations  are  required  to  clarify 
the  differences  in  t^ehavior  between  grades  of  graphite  having  high  ard  low 
specific  reactivity,  numerical  calculations  will  be  performed  for  "fast" 
and  "slow"  heterogeneous  reactions,  respectively.  In  particular,  in  the 
illustrative  examples,  we  will  utilize  eqs.  (4)  and  (b),  with  a  value  of 
n  =  1/2.  The  values  of  the  rate  data  will  be  arbitrarily  taken  as, 


(  E  =  44.  0  k  cal.  /mole 
i  k  =  6.  720  X  10*  lb.  /ft.  ^^^sec.  atm.^/‘ 


"slow' 


E  =  42.3  k  cal. /mole 
k  -  4.  47  3  X  10^  lb, /ft.  sec.  atm.^^  " 


In  the  rate  controlled  regime,  the  rate  of  oxygen  consumption 
adjacent  to  the  surface  is  so  low  that  tne  mass  fraction  of  the  element  oxygen 
will  be  essentially  the  saine  as  m  the  undisturbed  stream.  However,  the  atomic 
species  will  recombine  in  tne  low  gas  phase  temperature  regime  adjacent  to 
the  surface  and  hence  the  mole  fraction  of  molecular  oxygen  at  the  surface 
will  be  very  nearly  equal  to  its  value  in  undissociated  air.  Further,  the 
static  pressure  gradic”'  in  the  reaction  zone  =«  direction  normal  to  the 
surface,  is  negliu’ibte.  1  hese  considerations  therefore  enable  us  to  write. 


"O. 


V 


I'.r  iiif  rt-.iction  rat.-  coiitroUed  rt-sime  only,  and  hence  equation  (4)  becomes 


(*-.v)React.  ,  1/2  -E/RT 

kXo2 


1/2 


(W) 


where  Xqu  =  0-  21  for  undissociated  air. 

”w 

It  is  noted  if  oxygen  does  reach  the  surface  in  the  atomic  Jtate, 
reactions  will  occur,  but  at  a  somewhat  different  rate  (Ref.  57) ,  which  will 
promote  a  shift  in  the  transition  regime. 


The  cons.-rvatiLin  c()Udliun  for  species  L  bt  eomes: 


Tin-  I  oinpoiienl  i  •'  momentum  bei-omes: 


b)  Trt* ••.sport  and  Tbermodynamic  Properties 

The  transport  troefficient.s  requiretl  for  the  definition  of 
the  physical  problem  include  the  coefficient  of  viscosity  for  each  of 
the  n  pure  species,  (n  -r.)/2  symmetric  binary  diffusion  coefficieulo, 
and  n  self  diffusion  coefficients.  These  properties  may  be  calculated 
from  the  followini;  equations  (ref,  58  ). 


(24) 


/i 


S'-:  - 


5  -y/  IT  m  UT 


_3 

lb 


|T,,.3^3 


ITT 


(2b| 


v^h'.i  li  in-  till.-  r!;;<irf)ii-  i-Mu-li*  thvf>ry  I'oriniil.n-  for  the  vim  osity  ;intl 
tlu'  <iiffn-.u>i!  I  livfiy.  Iti  ortiur  to  evaluate  thost; 

urope  rti«->  on«:  re<iuir«;s  a  krjow  ie'l^v  ol  the  i  ollision  dianietur  0^, 
.iiiti  the  eolliMion  mtear.il  a  ■>  .i  funv  tion  of  the  reduced 
temp*" nature,  where  T^=KT/6. 

In  the  above,  the  symbol  li  '  is  the  reduced  mass  givuin  by: 


y.j 


'J 


<T''j 


and  the  coiitsion  diameter  by: 


<r,j  = 


(26) 


(27) 


'.'ote  further  that  the  constants  to  be  utilized  in  the  above  equations  are 
^iven  tn  Table  IV  and  the  resulting  transport  piuperlies  of  the  pure 
species  which  are  based  on  the  Rigid  Sphere  and  the  Lenrtard  Jones  6:12 
potential  respectively  are  shown  in  Figures  2,3  and  4. 
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c)  Similarity  Transformation 

Upon  introducing  the  Mangler-Oorodnitsyn  transformations 


^  ^>y^w^er^o(x  (29) 


and  assuming  that  local  similarity  holds,  then  equations  (19)  through 
(21)  may  be  reduced  to  a  set  of  ordinary  non-linear  equations.  The 


diffusion  equation  becomes; 


Is  -it  1^  =  0 

LVr  M  S  ^  f 


The  ccr.Jcrvation  of  momentum  becomes: 

'V] ' 


O  (31) 


while  thj  energy  equation  becomes: 


(fLi V  Cf  -f  S  [if  ^Jt] 


27 


d)  Chemical  Constraints  and  Boundary  Conditions 

It  will  be  assumed  here  that  the  total  number  of  dominant  Species 
present  in  the  nas  phase  is  six,  including  atomic  and  molecular  oxygen 
atomic  and  molecular  nitrogen  and  the  combustion  products,  carbon 
monoxide  and  carbon  dioxide.  Thus,  there  are  six  unknown  chemical 
source  terms  Wj  and  six  unknown  concentrations  Xj  at  each  point 
within  the  boundary  layer. 

At  the  surface,  there  are  six  unknown  concentrations  which 
must  satisfy  the  chemical  constraints  imposed  by  the  surface  oxidation 
process. 

Let  us  therefore  consider  the  chemical  constraints.  The 
conservation  of  r>>«»mical  species  in  the  gas  phase  requires  that 

vi  i  (33) 

i 

Since  it  has  been  assumed  that  the  nitrogen  species  do  not  react  with 
oxygen  or'  carbon  to  form  permanent  compounds  (i.  e.  the  formation 
of  NO  or  CN  is  not  considered  explicitly) ,  one  may  also  write  for  the 
case  of  thermochemical  equilibrium 
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while  combustion  stoichiometry  requires 


'^0  ^  '^COj  (35) 

^CO  ”  “  ^ 

For  the  case  of  frozen  reactions  one  has 

-SO  (37) 

The  surface  boundary  conditions  on  velocity  and  temperature  may 
then  be  written: 


The  boundary  conditions  on  the  composition  are  obtained  from 
the  simultaneous  solution  of  the  equilibrium  constants: 


2 


^0  Vi#  Xtow 

It  is  also  noted  that  in  the  theoretical  models  considered  here,  the 

nitrogen  atoms  are  permitted  to  recombine  on  the  surface. 

one  may 

also  write: 

)w  = 

and  hence  since 

=■ 

Lf  )-v 

(«) 

one  immediately  obtains  for  the  oxidation  rate 

1  - 

■  ^COv)v 

The  boundary  rnnditions  at  the  outer  edge  of  the  boundary  layer 


‘or  velocity  and  temperature  are  given  by; 


unti  in  addition,  one  has  n-l  relations  of  the  form: 

X ,  =  X 

The  total  number  of  boundary  conditions  therefore  equals  the 
order  ot  the  mathematical  system. 


V, 


DtsrfjSSIOX  OF  KFSIILTS 


In  this  section,  the  results  of  the  theoretical  investigation  will  be 
discussed  in  detail.  As  already  noted,  a  convenient  coordinate  system 
for' the  study  of  the  liypersonic  laminar  boundary  layer  is  tlic  bcdy- 
oriented  coordinate  system  shown  in  Figure  1.  In  order  to  treat  the 
problem  in  the  absence  of  experimental  data,  the  high  temperature  trans¬ 
port  coefficients  and  thermodynamic  properties  of  the  reacting  species 
Were  determined  theoretically  utilizing  kinetic  theory,  statistical  thermo¬ 
dynamics  and  the  gas  constants  given  in  Table iT.  The  theoretical  results 
are  shown  in  Figures  J  to  o.  Here  it  is  remarked  tnat  the  properties  of 
the  gas  mixtu'e  have  been  treated  as  variable  and  were  computed  through 
the  boundary  layer  as  a  function  of  both  the  local  gas  composition  and  the, 
properties  of  the  pure  species  at  the  local  gas  temperature.  Thus,  since 
the  gas  properties  are  computed  as  part  of  the  solution,  one  dots  not  require 
simplifying  assumptions  such  as  constant  Prandtl  and  Lewis  numbers,  or 
a  constant  product  of  density  and  viscosity  in  order  to  solve  the  problem. 
That  is,  once  an  appropriate  inlermolecalar  force  law  between  a  pair  of 
molecules  has  been  selected,  one  may  immediately  calculate  uniquely  the 
properties  the  pure  species,  and  then  one  may  proceed  to  utilize  these 
to  determine  the  mixture  properties  as  dictated  by  the  chemistry  of  the 
particular  physical  problem. 

In  Figure  7  are  shown  typical  values  of  the  variation  of  the  normalized 
product  of  density  and  viscosity  through  the  boundary  layer,  as  a  function  of 


the  stretched  norrrfal  cccirdinatc.  Qualitatively,  the  trends  are  the  same, 
namely  that  the  quantity  decreases  with  increasing  gas  tempera¬ 

ture.  This  may  readily  be  explained  as  follows.  In  a  dissociating  gas 
boundary  layer,  which  is  locally  a  constant  pressure  layer,  the  density 


ratio  is  given  by: 

.-A  =  .A  C4,^ 

A  iww 

and  hence,  clearly,  the  decrease  in  molecular  weight  with  increasing  temp¬ 
erature  produces  a  smaller  density  ratio  (Sweater  density  change) 

with  increasing  gas  temperature  than  will  occur  in  a  non-dissociating  gas. 
This  .8  particularly  true  for  the  case  of  the  mass  transfer  of  foreign 
species  at  the  surface  having  a  higher  molecular  weight  than  the  primary 
gas  species.  Since  the  increase  in  viscosity  witn  gas  temperature  is  less 
than  linear,  the  increase  in  the  viscosity  ratio  with  increasing 

temperature  does  not  compensate  for  the  decrease  in  the  density  ratio 

•  Thus,  the  larger  the  ratio  of  the  gas  temperature  at  the  edge  of 
the  boundary  layer  to  that  at  the  wall,  the  smaller  will  be  the  value  of  the 
normalized  product  at  the  outer  edge  of  the  boundary  layer,  and  the 

poorer  the  assumption  that  constant.  This  is  shown  in  curves  1, 

2,  and  5  of  Figure  7,  which  were  calculated  for  the  chemical  constraint  of 
loC'.l  equilibrium  flow.  By  way  of  further  comparison,  curve  3  of  Figure  7 
was  computed  for  the  case  of  equilibrium  dissociated '.air  (no  graphite  com- 
buscion) ,  and  curve  4  was  calculated  for  tne  chemical  constraint  of  frozen 


reactions,  all  for  the  same  flight  conditions.  It  in  seen  that  under  certain 
conditions  the  density-viscosity  variation  may  have  nearly  the  same 
behavior  for  both  the  frozen  flow  and  local  equiUbriunrt  constr;>ints, 

(  compare  curves  4  and  3) ,  although  the  other  gas  properties:  surh  as,  for 
i!xample_  the  thermal  conductivity,  do  not. 

When  one  defines  the  F^randtl  number  of  the  gas  in  terms  of  the 
ti-o»en  specific'heat  of  the  mixture  Cp  ,  the  viscosity  of  the  mixture yU 
and  the  frozen  thermal  conductivity  K,  one  obtains  the  results  shown  in 
Figure  H.  Here,  the  value  of  the  Prandtl  inimiier  of  the  gas  at  the  surface 
versas  the  surface  temperature  over  a  aide  range  of  pressures.  It  appears 
from  this  figure  that  the  presence  of  carbon  dioxide  at  the  surface  at  the 
lower  surface  temperatures  acts  to  increase  the  Prandtl  number.  At  a 
given  value  of  the  temperature,  the  effect  of  a  pressure  increase  is  to 
decrease  ihe  dissociation  of  carbon  dioxide  and  hence  the  shifting:  equil¬ 
ibrium  composition  with  increasing  pressure  results  in  an  increase  in  the 
Hrandtl  number.  The  total  effect,  however,  is  less  than  four  percent 
si.ice  th  r-  x.'mum  mole  fraction  of  carbon  dioxide  is  less  than  0.  i.  Of 
grcd*.  r  interest  are  the  Lewis  numbers  which  have  hee.a  defined  here  in 

terms  of  the  multicomponent  difiusion  coefficients  D..,  which,  unlike  the 

tl 

binary  difiusion  coofficionts  •  dire  not  Rvmuxtric.  Consequently, 

since  all  other  gas  properties  which  appear  in  the  multicomponent  Lewis 
numhiir  .•■re  ntit  dependent  on  tt.e  particular  p.air  of  particles  involved  in  a 
given  binary  enedimter,  while  the  multicomponent  diffusion  coefficients 
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arc  both  temperature  and  composition  dependent.  It  is  clear  that  there 
2 

are  cenerally  n  -n  values  of  the  Lewis  number;  in  this  case  there  are  30, 
of  which  12  representative  values  are  shown  in  Figure  9. 

When  the  boundary  layer  equations  arc  integrated  subject  to  the  appropriate 
hoiMidary  conditions,  one  obtains  the  velocity  distribution,  the  temperature 
distribution  and  the  distribution  of  species  through  the  layer. 

These  are  shown  in  Figure  10  and  11.  Corresponding  chemical  source 
terms  are  shown  in  Figure  12. 

Since  preferential  diffusion  can  occur,  it  is  found  that  the  ratio  of  the  element 
oxygen  to  the  element  nitrogen  at  the  surface  is  not  necessarily  equal  to  its 
value  at  the  edge  of  the  boundary  layer,  and  is  in  fact  an  eigenvalue  of  the 
mathematical  system.  This  is  shown  in  Figure  13  and  14. 

In  the  rate  controlled  regime,  it  is  found  that  the  ratio  of  the  mass  fraction 
of  cai  uun  monoxide  to  the  mass  fraction  of  carbon  dioxide  is  a  sensitive  function 
of  surface  temperature,  which  increases  with  increasing  surface  temperature. 
This  can  be  shown  to  be  the  case  from  both  theoretical  and  experimental  con¬ 
siderations  and  is  demonstrated  in  Figure  15. 

Examination  of  the  literature  indicates  that  the  reaction  rat,,  controlled 
moss  transfer  is  also  a  very  sensitive  function  of  surface  temperature  and  in 
fact  as  many  different  activation  energies  and  specitic  reactivities  have  been 
reported  in  the  literature  as  there  are  combinations  of  investigators  and  grades 
of  graphite. 
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The  work  of  five  different  investigators  is  shown  in  Figure  16.,  The 
e\[)<jnenlial  dependence  of  the  oxidation  rate  on  surface  temperature  is  clearly 
si-en  which  leads  to  the  usual  Arrhenius  representation. 

Typical  values  of  ihe  reaction  rate  results  therefore  appear  in  Figure  17 
for  «  representative  flight  condition. 

In  Figure  '8  we  have  shown  the  correlated  results  for  the  mass  transfer 
of  oxidation  products  for  both  the  rate  controlled  and  diffusion  controlled 
regimes.  When  the  surface  temperature  exceeds  J000°R  the  oxidation  rate 
levels  off  and  becomes  insensitive  to  the  magnitude  of  surface  temperature, 
simply  because  the  mass  loss  is  controlled  by  the  diffusion  of  oxygen-bearing 
speci  -s  to  the  surface  rather  than  the  specific  reactivity  of  graphite. 

Upon  making  use  of  the  concept  of  th*-  sum  of  resistances  in  scries, 
one  m<iy  write: 


I 


n-.w 


React. 


(bU) 


where  the  reaction  term  predominates  at  low  surface  temperatures  and  diffu¬ 
sion  predominates  at  high  surface  temperatures.  This  equation  was  utilized 
in  calculating  the  transition  regime  in  Figure  19. 

In  order  to  obtain  an  over-all  feeling  for  the  combined  effect  of  pressure  and 


surface  temperature  on  the  oxidation  rate,  Figure  ZO  has  been  prepared.  Here 
it  is  seen  that  th're  are  4  different  regimes  for  the  oxidation  of  graphite.  In 


rate  was  correlated  by  means  of  the  foUowinf{  equation: 


Pe  f(A) 


BTU 

=  33.  3  .  0.0333(H,.h„^^).  . . 

ft.  sec.  atm. 


I  where  ;(A)  represe:as  me  effects  of  yaw,  ^  is  the  kroneclcer  q  such  'hat 

thi<?  equation  applies  to  both  axially  symmetric  (  ^  =  1)>  and  two-dimensional 

(  S~  geometries.  P^,  is  the  local  static  pressure  at  the  outer  edge  of  the 

boundary  layer,  H  is  the  stagnation  enthalpy  and  h  is  the  static  enthalpy 
*'  Air  . 

of  the  gas  evaluated  at  the  surface  temperature  and  treated  as  if  it  were  air, 
expressed  in  BTU /lb.  Figures  23  and  24  represent  different  ways  of  depicting 
the  correlated  heat  transfer. 

If  one  utilizes  the  horizontal  portion  of  the  graph  shown  in  Figure  25 
for  m^,  one  obtains  the  following  result  for  the  diffusion  controlled  mass  loss: 


’w  Diff. 


=  6.  2  X  10' 


3/2  1/2 

ft.  sec.  stm. 


Note  that  equation  (14),  which  represents  the  reaction  rate  controlled 
mass  loss  is  independent  of  the  geometric  factors  whereas  the  diffusion  con¬ 
trolled  mass  loss,  equation  (S3),  i.not. 

It  is  noted  that  in  both  equations  (S2)  and  fSl),  the  effect.s  ut  yaw,  which 
act  to  decrease  the  heat  and  mass  transfer,  appear  in  the  function  f(A).  Although 
more  cornplicated  correlations  have  been  suggested,  (Ref.  S9),  the  author 


believes  that  the  following  approximation  is  s.atisfactory. 
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It  is  of  some  interest  to  calculate  the  so-called  "effective  heat  of 
ablation".  This  may  be  done  readily  for  the  diffusion  controlled  regime  by 
taking  the  ratio  of  to  ,  and  one  obtains  a  result  which  is  independent 
of  the  geometry  of  the  body,  that  is: 

Q*  5  570  +  5.  37  -  h^^  ).  BTU/lb.  (54) 

•  li  r 

It  should  be  noted  that  since  the  driving  force  for  the  mass  transfer 
process  is  nrit  really  the  heat  transfer  but  rather  the  oxidation  process,  one 
should  not  be  surprised  that  the  effective  heat  of  ablation  is  of  the  order  of 
70,000  BTU/lb.  at  orbital  velocity.  It  must  be  understood  that  graphite  is 
actually  an  excellent  sink  material  which  also  undergoes  surface  oxidation. 

In  this  paper  we  have  also  determined  the  surface  shear  stress 
defined  by 

Tw  =  /w  <'It>w 

and  a  typical  result  is  shown  in  Figure  26. 

Upon  defining  a  local  skin  friction  coefficient  c,  : 


and  a  Reynolds  number  based  on  x: 


one  also  obtains  the  correlated  results  shown  in  Figure  27.  For  purposes 
ot  comparison,  results  arc  also  given  for  the  skin  friction  coefficient  in  the 
absence  of  surface  oxidation. 

.  Finally,  it  may  be  of  some  interest  to  be  aole  to  rapidly  estimate 
the  surface  temperature  as  a  function  of  environmental  conditions  such  as 
Hj,  and  P^.  In  general,  one  has: 


where  is  the  total  convective  heat  flux  to  the  surface,  (T*  €  represents 

radiaiio.t  to  li.v  surface  from  the  hot  gas  cap,  (P **  reradiation  from 
“O  T 

the  surface  and  (K  )  ,  is  the  heat  conducted  into  the  interior  of  the 


3  y  solid 


solid  graphite. 


If  one  assumes  that  radiation  equilibrium  obtains,  the  c'^ht  hand 
side  of  equation  (98)  may  be  set  equal  to  zero.  In  this  case,  one  may  write 
that  the  radiation  equilibrium  temperature  is: 


r  " 

\  OTi 


<re„T 


In  general,  for  suborbital  flight,  0"tf  T  ^<.  Q  and  if  this  term  is  neglected 

S  g  w 

one  obtains  the  typical  results  shown  in  Figure  28  for  a  value  of  the  radiation 

parameter  ^  J  R  '  =  0.  85  ft. 
w  *1  B 


41 


VI.  CONCLUSIONS 


A  study  of  the  ablation  of  graphite  indicates  that  one  may  recognize 
at  least  four  distinct  acrothermochemical  regimes.  As  the.  surface  temper¬ 
ature  rises  into  the  threshold  range,  the  mass  transfer  process  is  initially 
rate  rnntrclled  ^nri  ollows  a  fractional  (half  order)  reaction  law.  At 
surface  temperaturec  of  approximately  1400°  R  <  <  3200°  R,  a  transition 

regime  exists  and  the  overall  rate  of  the  process  is  controlled  both  by 
cherrical  and  gas  dynaniic  factors.  .4t  higher  surface  temperatures  extend¬ 
ing  til  rough  6000°  R,  the  process  is  diffusion  controlled  and  the  rate  of  mass 
loss  from  the  surface  is  limited  by  the  rate  at  which  oxidizing  species 
diffuse  to  the  surface.  Finally,  there  is  a  regime  in  which  the  surface 
temperature  ia  sufficiently  high  so  that  the  sublimation  rate  exceeds  the 
surface  oxidation  rate,  and  carbon  atoms  and  molecules  are  present  in  the 
gas  phase. 

In  the  rate  controlled  and  transition  regimes,  the  ablation  rate  ot 
graphite  is  a  sensitive  function  of  the  surface  temperature,  and  depends 
critically  on  the  specific  reaction  rate  which  has  an  exponential  temper¬ 
ature  dependence.  It  is  in  the  low  surface  temperature  regime  that  one 
can  distinguish  between  different  grades  of  graphite,  since  the  activation 
energy  and  the  effective  collision  frequency  depend  on  the  molecular 
structure  of  the  particular  carbonaceous  material.  Therefore,  unless 
one  has  an  independent  knowledge  of  the  chemical  kinetics  of  the  combustion 
process  for  the  specific  carbonaceous  material,  one  cannot  make  a  precise 

■l£ 
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prediction  of  the  oxidation  rate  for  surface  temperatures  below  3200°  R. 

1  ne  thermocheniical  response  of  graphite  in  the  transition  regime 
can  he  easily  synthesized  from  a  knowledge  of  its  behavior  in  the  purely 
rate  controlled  and  diffusion  controlled  regimes.  A  more  precise  study 
of  the  transition  regime  itself  does  not  appear  to  be  warr:."teil  at  this  time 
beca  ise  of  an  incomplete  knowledge  of  certain  of  the  chemical  factors,  e.  g. 
the  CO/CO^  ratio.' 

In  the  diffusion  controlled  regime,  it  is  found  that  the  ablation 
rate  is  proportional  to  t.-ie  square  root  of  the  pressure  and  is  essentially 
independent  of  surface  temperature.  It  is  further  noted  tiiat  for  surface 
temperat  ires  above  3200° R,  barring  spalling  effects,  the  mass  loss  Is 
theoretically  independent  of  the  grade  of  the  graphite.  The  ma  ‘f  uncertainty 
resides  in  the  somewhat  arbitrary  selection  of  the  high  temperature  transport 
properties  due  to  the  lack  of  experimental  data  for  any  of  the  six  gaseous 
species  present  at  the  temperatures  of  interest.  However,  judging  by  a 
comparison  between  theory  and  experiment  for  stagnation  point  heat  transfer 
in  the  absence  of  mass  transfer,  it  is  estimated  that  the  uncertainty  in  the 
transport  properties  introduces  no  more  tlian  approximately  a  fifteen  percent 
error  in  tne  present  final  results. 

A  study  of  the  calculated  results  indicates  that  in  the  diffusion 
controlled  regime,  the  h-sat  and  mass  transfer  are  independent  of  the  rates 
of  the  homogeneous  reactions,  provided  that  the  heterogeneous  reactions  are 
rapid.  That  is,  when  it  was  assumed  that  thermochemical  equilibrium  was 


M 


obtained  at  the  surface,  but  not  necessarily  throut{hout  the  boundary  layer, 
bot^*  the  "frozen  reactions"  and  "local  equilibrium"  chemical  constraints 
yielded  esi^eniially  identical  results. 

.  In  considering  the  heat  conducted  into  the  solid,  one  notes  that 
the  separate  contrU-  ri&ns  of  convection,  diffusion,  conduction  and  radiation 
must  be  included.  In  this  study,  it  was  found  that  the  net  heat  transfer, 
including  combustion  effects, exceeds  tne  aerodynamic  heat  transfer.  One 
may  conclude  therefore,  that  the  exothermicity  of  the  oxidation  reactions 
exceeds  the  decrease  in  heat  transfer  due  to  mass  transfer  cooling  effects, 
(thickening  of  the  boundary  layer),  which  results  in  a  net  increase  ir.  the 
heat  transfer  to  the  solid.  However,  this  increase  is  usually  less  than 
ten  percent  for  the  full  range  of  interest. 

With  regard  to  viscous  drag  effects,  it  is  concluded  that  for  the 
hypersonic  flight  regime  the  reduction  in  skin  friction  due  to  the  mass 
transfer  of  wxidation  products  does  not  appear  to  be  BigiiiTicant,  and  is  of 
the  order  of  seven  percent  or  leee. 
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